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We investigate the survival probability of a localized 1-d quantum particle subjected to a time 
dependent potential of the form rU (x) sin cot with U (x) = 28(x ~ a) or U (x) = 2S(x — a) — 2S(x + a). 
The particle is initially in a bound state produced by the binding potential —2S(x). We prove 
that this probability goes to zero as t — » 00 for almost all values of r, cj, and a. The decay is 
initially exponential followed by a t~ 3 law if u is not close to resonances and r is small; otherwise 
the exponential disappears and Fermi's golden rule fails. For exceptional sets of parameters r,uj 
and a the survival probability never decays to zero, corresponding to the Floquet operator having 
a bound state. We show similar behavior even in the absence of a binding potential: permitting a 
free particle to be trapped by harmonically oscillating delta function potential. 



PACS: 03.65.Db; 03.65.Ge; 32.80.Fb 

1. Introduction 

Quantum systems subjected to strong external time 
dependent fields often show very complex behavior, e.g. 
the ionization probability of an atom may be a compli- 
cated function of the frequency, amplitude, pulse shape 
and other parameters of the field [1-4]. Such phenomena, 
which go beyond conventional perturbation theory in the 
field strength (r here), are readily seen in numerical solu- 
tions of the time dependent Schrodinger equation. There 
are also various approximate analytic methods which re- 
produce many experimental results [5,6,7], but there is 
however no rigorous theory of such phenomena even for 
model systems. In this paper we describe new exact re- 
sults for a toy model which has both bound and contin- 
uum eigenstates, subjected to a harmonically oscillating 
potential. They reveal a very rich structure for the time 
evolution of even very simple quantum system. In partic- 
ular the transition from a bound state to the continuum 
is seen to be much more complex than the simple expo- 
nential decay obtained from conventional perturbation 
theory via Fermi's golden rule [2]. 

To obtain exact results we need to consider simplified 
model systems. In particular we cannot treat (for the 
present time) realistic description of the interaction be- 
tween radiation and matter, such as the dipole approx- 
imation. We note however, that a comparison of our 
earlier results on an even simpler version of this model 
with experiments on the strong field ionization of Ryd- 
berg atoms [8] showed some surprising similarity between 
the two. We interpret this as an indication of a certain 
universality in the Schrodinger evolution of a system with 
bound and continuum spectrum subjected to time depen- 
dent external forces and that our model retains some of 
its behavior. 

2. The model 

The important idealizations in our model are: 1) Space 



is one-dimensional, 2) The "internal" potential creating 
the bound state is given by an attractive delta function 
at the origin, 3) The interaction with the external field 
has the form rj(t)U(x) where r) is periodic in time with 
a rectangular envelope, and 4) U(x) is given by one or 
more delta functions at different locations on the x-axis. 
The first two assumptions are quite common for model- 
ing short range binding potentials [9,10] and should not 
affect greatly the basic physics of the ionization process. 
Assumption 3 means that we do not consider situations 
[11] where there is some "ramping" in turning the exter- 
nal field on and off. This should not be too serious when 
the pulse is of a long duration compared to the period of 
the field which is the case we are concerned with here. 
Assumption 4 on the other hand clearly makes the inter- 
action in our model very different from the real interac- 
tions between radiation and matter: dissociation due to 
electromagnetic fields are described approximately by a 
dipole interaction of the form U(x) — x. Unfortunately 
we have not been able to obtain exact results for this 
case beyond those described in [11]. The only feature of 
the dipole interaction we are able to mimic is the spatial 
symmetry. 

Using suitable units in which H = 2m = 1 (m is the 
particle mass) the time evolution of our system is given 
by the Schrodinger equation 



dx 2 



28(x) + U(x)r)(t) 



iP(x,t), (1) 



where r/(t) = rsinwt and the parameters r,oj represent 
the amplitude and frequency of the time dependent po- 
tential. The spatial structure of the external potential 
will be taken in two forms: U%(x) = 26 (x — a) and 
U2^x) = 2[S(x + a) — 8(x — a)]; U2 has the symmetry of the 
dipole interaction. The factor 2 in front of the binding 
potential is chosen so that the unique bound state of the 
unperturbed system is Ub(x) = e~' x ' with binding energy 
Eq = luq = 1. Eq.(l) is to be solved subject to the initial 
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condition ip(x,0) — Ub(x). We can readily extend our 
methods to more general sums of delta functions. Our 
main interest is in the survival probability of the bound 
state at time t: \9{t)\ 2 = \(tp(x,t)u b (x))\ 2 . 

The case U\(x) with a = 0, which corresponds the 
parametric perturbation of the binding potential, was 
treated in [12,13]. We showed there that |#(i)| 2 has 
both exponential and power law parts which are well 
separated only when the strength of perturbation r is 
small. This was true for all uj away from resonances 
(uj ,6 A^ -1 , N an integer), with uj < 1 corresponding 
to ionization via "multi-photon" processes. We also ob- 
tained there non-monotonic dependence of the escape 
rate on r and uj. Somewhat to our surprise we found 
qualitative (and even semi-quantitative) agreement be- 
tween the predictions of this model about resonance 
behavior of the survival probability of localization and 
some experimental observations on the ionization of Ry- 
dberg atoms by strong microwave fields [8]. It was also 
proved in [13] that when r](t) is a sum of a finite number 
of harmonics, {^2j—i Aje l ^ ut + complex conjugate), then 
the survival probability |#(t)| 2 — > as t — > 00 for any 
M < 00, Am 7^ 0. There are however very special infi- 
nite sequences Aj, given explicitly in [13], for which we 
proved that the system never ionizes fully |#(i)| 2 ~h 0. 

Here we show that the situation is quite different and 
much richer when we consider U\(x) with a / or 
U2(x). (In physical units the position of perturbation 
corresponds to ha/\/2rnEo). In particular we prove that 
for rj(t) = rsinwi there exist two-dimensional manifolds 
in the space of the three parameters uj, r, a on which 
|#(t)| 2 ■/* as t —* 00. (We shall take without loss of gen- 
erality uj, r, a positive.) This means that while |#(t)| 2 — > 
for almost all parameter values of the forcing, "excep- 
tional" cases can also be constructed quite readily. This 
does not occur for fixed uj and a if r is small enough and is 
thus outside conventional perturbation theory. We find 
in addition that when uj is very close to a resonance, 
uj w A r ~ 1 + dynamic Stark shift, then the decay may not 
have the exponential part predicted by the golden rule 
no matter how small r is even when |#(i)| 2 — > 0. 

We can also consider the case when there is no binding 
potential at all, i.e. the term — 26 (x) is absent in (1). 
In this case we have for rj(t) = a free particle, which 
when initially localized in the vicinity of the origin will 
diffuse away: the probability of being in any fixed region 
decaying as i -1 . On the other hand the perturbations 
with special values of a,uj,r can make the particle stay 
localized for all time. 

3. Results for U 2 (x) 

We give here an outline of the proof which follows along 
the lines presented in detail for the case a = in [13]. 
Expanding ip(x,t) in terms of the eigenf unctions of the 



unperturbed Hamiltonian H = — ^ — 26 (x), we write 

POO 

ip(x,t) =6(t)e-W +it + / e(k,t)u(k,x)e- lk2t dk, (2) 



where the initial conditions are 8(0) = 1, O(fc,0) = 
and the explicit expression [12] for the continuum states 
are 



u(k,x) 



1 



Akx 



i\kx\ 



i + i\k\ 



-00 < k < 00. 



Substituting (2) into (1) we obtain tp(x,t) in the form 
of a functional of tp(a,t) and ip(—a,t) while 6(t) can be 
written as 



9{t) 



l + 2ie~ a [ ri(s)[Y + (s)-Y-(s)]ds, 
Jo 



(3) 



where Y ± (t) 



t(j(±a,t). We now express <d(k,t) in 



terms of Y ± (t) too, take x — ±0 in (2) and obtain a 
coupled pair of integral equations 

Y ± (t)^e- a +f 7 1 {s)[K ± it^s)Y+{s)-K^it- s)Y~(s)] ds. 
Jo 

(4) 

The Laplace transforms of the kernels K ± (t) are 



k-(p) = 



VI - ip - 1 



a/1 - ip 



with the choice of branch a/1 — ip — * 1 when p — > 0. For 
rj(t) = rs'mujt, letting y ± (p) be the Laplace transforms 
of Y ± (t) and setting f(p) = y + (p) — y~{p), eq.(3) yields 



9(0 



2-ki J c p 



iuj) - f(p + iw)}dp (5) 



with integration along a contour C which goes from —zoo 
to too in the right half plane avoiding p — by a small 
semi-circle in the left half plane. 

Survival of bound state: The survival probability \6(t) | 2 
is determined (see (5)) by the analytic structure of y ± (p). 
Setting y^ = y(p + iuin), k^ — k^ip + iujn) eqs.(4) 
turn into the recurrence relations for the vectors y n = 



Vn = 



ir 
~2 



k+ 



(jM-l-JM+l), (6) 



p + turn 

where y n , fc* may be viewed as functions of the param- 
eter p in the strip $s(p) € [0, uj). The poles of y(p) are in 
the left half-plane at £ + ium, £0 < 0, and the branch 
points at p = — i — inuj, n € Z (the latter ones are in- 
herited from k^(pj). After making horizontal cuts at 
p = x + inuj, —00 < s = SRp < we push the contour C 
in (5) along the branch cuts into the left half plane. The 
poles then contribute a series of residues with the com- 
mon exponential factor e^°* while the integrals around 
the cuts generate a contribution in the form of a series in 
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terms of t^^ 1 / 2 , j > 1 (see [12]). The imaginary part 
of £o> 3£oi is identified as the dynamical Stark shift [2,7] 
of the resonance frequency and T — — 23?^o is the decay 
exponent in the initial stage of evolution when r is small 
and the exponential and polynomial parts of |#(i)| 2 may 
be separated. 

If 3?^o = 0, i.e. the poles lie on the imaginary axis, 
then 9{t) -f* as t — > oo. This can happen when on the 
imaginary p axis the homogeneous recurrence, 



(z„- 



Zn+l) 



(7) 



fc+ k n 

2 — , 

associated with (6), has non-trivial solutions which decay 
sufficiently rapidly as n — > ±oo, i.e. l z n| 2 < oo. This 
is a manifestation of the Fredholm alternative [14]. We 
show now that unlike the case a = treated in [12,13] 
such solutions though non-generic are possible for U 2 (x) 
and for U\{x) with a > 1/2. 

Setting p — ig (with a real g) we construct a particular 
solution of (7) for w > 1. It is clear that if Zj = Zj+i = 
then all successive z n will be zero too until the matrix in 
(7) becomes degenerate. We set z n = for all n < and 
require the determinant (k^) 2 — (k^) 2 to vanish, which 
allows z\ 7^ 0, in particular z\ = zf . This implies 

ay/-l-go = itN, N = 1, 2, (8) 

where the parameter go e (— ui, —1) represents the "bind- 
ing energy" of a localized state produced by the pertur- 
bation (Floquet state [6,7]). For n > 1 all k„ are real 
positive and the matrices in (7) are non-degenerate. By 
inverting them it is easy to show that z~ = (— l) ,l+1 z+ 
and obtain a scalar recurrence for z+, n > 0. Using a 
new variable p n = iz^ +1 /z+ the recurrence (7) takes the 
form 

2 1 2 

Pn = —j , Pi = ~r, (9) 

rk n p n _! rki 
where n > 1 and fc„ = £;+ + (— l) n+1 fc~. 

A careful analysis for the case N = 1 in (8) shows that 
the decaying solutions of (9) can be constructed with a 
unique r = r(a,co). Two inequalities, 

kik 2 >2k 3 (2k 3 -k 2 ), hk 2 > fc 3 (4fc 3 -k 2 ), (10) 

which are necessary and sufficient for the existence of 
solutions of (9), specify regions in the a, u>, plane. Rela- 
tions (10) can always be satisfied if we choose go +ui <c 1 
which makes k\ large. The second of them gives the up- 
per bound on u> which becomes very strict, < Lu+go <C 1 
when —go ~ 1, i.e. a is large, see (8). The interval where 
the stabilizing r is located can be specified too and we can 
prove that for an arbitrary frequency ui > 1 there is an in- 
terval of a = 7r/ y/— 1 — go with go in (— u, —1) where, for 
a particular choice of amplitude r = r(a,u) the system 
does not ionize completely. Instead 0(t) — > e l9ot F(u> s t), 
where co s is the stabilizing frequency of the perturbation 
and F is a periodic function with period 2ir. 



4. Results for Ui(x) 



This type of stabilization of the bound state takes place 
for the perturbation with the potential U\{x) too, but 
only when a > 1/2. In this case there will again be a 2-d 
manifold in a,r,uj variables for which the bound state is 
stabilized and 9{t) — > e l9ot F(u s t), a quasi periodic func- 
tion of t. We also computed 6{t) numerically for this 
model by solving the integral equation for Y(t) and the 
most representative curves are shown in Fig.l. The slow- 
est decay on the time interval < t < 1607r/wo w 500 
was obtained for u) = 1.12 which is close to the value of 
uj s Rj 1.089..., evaluated by constructing the decaying so- 
lution of (9), and pushing it to as large n as we can within 
the accepted precision. |0(£)| 2 near the Stark shifted res- 
onance, ui = 1.2, has no interval in which the decay is 
exponential in contrast with such a decay observed for 
uj = 1.25 and uj = 0.8. 
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FIG. 1. Plot of log 10 



The ripples on the curves in Fig.l have the frequency 
of the perturbation lo and the modulation is due to beats 
among u>, u>o an d go- 

A perhaps simplistic explanation of the stabilization is 
maybe to view it as some kind of bouncing of the trapped 
particle between the delta potential wells. Note that the 
eigenfrequency go of the new bound state is a function 
(see (8)) of a only, but the amplitude r or the frequency 
ui must be fine tuned as functions of go to prevent the 
leaking out of the particle wave function. 

5. Time decay of the bound state 

For the model with Ui(x) the recurrence (6) has a 
scalar form and its solution can be written in terms of 
continued fractions which converge quite rapidly when r 
is small. For w > 1, neglecting terms of order of r 4 and 
higher, one may truncate the recurrence around each n 
by taking y n +m = if \m\ > 2. The solution of the 
truncated system for n = 0, which gives the main contri- 
bution to 6{t) has a pole at £o = 0(r 2 ), that solves the 
equation 

l + ^k+(to)[k+(Zo) + kl 1 (Zo)]=0(r 4 ). (11) 

Using (11) the contour integration in (5), where Y~ = 0, 
yields 
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6(t) 



g£o(w)t_|__ 



^[{uo 2 - (1 + A) 2 ][(w - 1 - A)t + l]3/2 ' 

(12) 

where Sft&M = -r 2 A(w)Vw - 1 - A, = A = 

r 2 cr(u;), and A(w),cr(w) 7^ are of order e~ 2a when a is 
large. A represents the Stark shift [6,7]. 

The survival probability |6*(t)| 2 has initially an expo- 
nential regime where it decays as e~ rt , with T = — 2K£o 
proportional to r 2 in agreement with Fermi's golden rule. 
As t increases an increasingly important role is played by 
transitions back to the bound state with probability pro- 
portional to the density near the origin as given by the 
second term of (12). As a result, for t S> the survival 
probability follows the power law decay |#(t)| 2 ~ i~ 3 , (see 
also [4,13,15,16]). The mathematical origin of this power 
law is the square root branch point at the bottom of the 
continuous spectrum. Note that this is much faster than 
when an initially localized free particle is permitted to 
evolve. The probability of it remaining localized then 
decays as t^ 1 [17]. 
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FIG. 2. log 10 |#(£)| 2 for two regimes of decay (plot of 
Eqs.(12),(14), note the different time scales and drastically 
different rates of decay). The upper curve for resonant ui 
does not follow Fermi's rule. The lower curve with the reg- 
ular behavior (exponential with a t~ 3 tail) is farther from 
resonance. 

6. Resonances 

Near the one photon resonance, u > 1 (11) implies 



£oM + r 2 A(l) v /^-l-e M = ir 2 a(l) + 0(r 4 ). (13) 

The solution of (13) gives A = 0{r 2 ) and T = o(r 2 ) 
whose dependence on r is determined by the order in r 
of u — 1, but in the case when w is very close to 1 + 
A we cannot separate anymore the contributions in (3) 
because the poles of y(p) are too close to the branch 
points. Integration in (5) yields 



0(t) 



iet+i-n/i roc g-tr 4 A 2 (l)x 2 74^,2^ 



(14) 



, [x 2 - h{a)} 2 + ix 2 ' 

where h{a) = 0(1) when a is not large. For r 4 t 3> 1 
cq. (14) implies |#(t)| 2 — > t~ 3 , but the initial exponential 



regime does not exist even for small r, see curves in Fig. 2, 
which are constructed using (12)-(14). (Compare with 
[4,15,16]). 

"Multiphoton" resonances. Let us compute the decay 
exponent when to < 1 which corresponds to a "multi- 
photon ionization" [6,8,12] for our simple model. To lo- 
cate the singular point p (which should be somewhere in 
the left half-plane very close to the imaginary axis) we 
require the homogeneous recurrence (a scalar analogue of 
(7)) to be solvable by a sequence z n such that \z n \ — > 
when \n\ — > oo. Using (9) and the continued fractions 
we represent po in terms of p^ 1 for n < and p n when 
n > 0. Both representations are rapidly convergent for 
r«l and the condition of their matching yields explic- 
itly the solution for £ in the lowest order A = 0(r 2 ) as 
before and T = 0(r 2N ). Zeros of fc_„ change the order 
of the decay exponent Y of the multiphoton ionization 
from r 2N to r 2N+2 which drastically slows down the the 
ionization rate. The stabilization can also be expected 
for lu w 1/2, 1/3, ... and therefore when u) < 1/2 we prac- 
tically cannot see the decay if t is not extremely large. 

7. Trapping of free particle 

Let us remove in (1) the binding potential and take the 
wave function 



/oo />oo 
F(k)u{k,x)dk, / \F(k)\ 2 dk = 
-OO J — OO 



1. 



which describes a localized particle in the vicinity of the 
origin. Here u(k,x) — (2Tr)~ 1 / 2 e lkx are the free eigen- 
functions. For r)(t) = the initial state evolves as 



i>\x,t) 



/CO 
F(k)u(k, x)e- lkH dk -» r 1 ' 2 , (15) 
-oo 



when t — > oo. Using in (1) the expansion in terms of 
u(k,x) and Laplace transform (ip(x, t) — > ip(x,p)) we ob- 
tain the infinite set of equations 



i> n (x) = i>° n (x)+ 

-\K{ x ){Vn-l - Vn+l) - k n( x )(Vn-l ~ Vn+l)l 



(16) 



where F n (x) means F(x,p + iu;n), = ip(±a,p + urn), 
and 



e i\x^a\^/ip 



ip — > i as p — > i. 



(17) 



Setting x = a and x = —a in (16) we obtain a recur- 
rence similar to (6) with vectors {i/>° (o), a )}- If on 
the imaginary axis the homogeneous recurrence (7) has 
a properly decaying solution for some a,u>,r then all y n 
have poles at p = ig(a) + icon respectively and therefore 
their inverse Laplace transforms which represent il>(±a, t) 
do not vanish as t — > oo. For a given lu the requirement 
that the determinant (fed) 2 — (fc^) 2 vanishes gives the pa- 
rameter g e (— co,0) as a function of a: g = —{■nN/a) 2 
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(N is an integer and clearly a > w/y/HJ). Setting z n = 
for all n < and z\ = z± , we can again invert the 
matrices in (7), obtain (9) and repeat the previous com- 
putation. The explicit form of the coefficients, 



I _|_ f -nn+lg-2aV9o+wra 



(18) 



implies here too a possibility to find a decaying sequence 
p n and a rapidly decreasing set of <~ r n u>~ n / 2 /Vn\, 
as n — ► oo. 

The poles of y„ develop by (17) into poles of ip{x,p) 
at the same points p = igo + ium. Therefore as t — > oo 
the wave function il>(x,t) will survive near the origin an 
can be represented as a series related to poles of (16) in 
Laplace space 



oo 

e igot ^ \e-\' x - a ^ 9a+un Q n {Ljt) + e-\ x+a ^ 9a+un P n (u)t) 

n=l 



where the coefficients P n , Q n are periodic functions which 
decay rapidly as n tends to oo. 

8. Concluding remarks 

Our results show the richness of the structure exhib- 
ited by a simple toy model driven externally in the pres- 
ence of a continuum. The survival probability can change 
greatly, including trapping in a localized state, as the pa- 
rameters of the external forcing arc varied. While not all 
the features of this simple model can be expected to be 
mirrored by real atoms driven by electromagnetic fields 
we believe that some features are rather universal. In 
particular the power law decay [16] and the Fermi golden 
rule violation at resonances are expected to occur quite 
generally. Even if the location of resonance is not on the 
real axis of the energy plane but has a small imaginary 
component and therefore the localized state is slowly de- 
caying, the power law tail can compete with the exponent 
on the whole observable time interval. 
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